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FOREWORD 

The  area  of  Sonobuoy  performance  prediction  is  of  vital  Importance  for  continued  im¬ 
provement  in  our  Antisubmarine  Warfare  (ASW)  capabilities.  As  part  of  the  Naval  Ocean¬ 
ographic  Office's  investment  in  ASW.  the  Environmental  Systems  Division  specializes  in 
on-scene  predictions.  This  report  identifies  the  major  on-scene  sonobuoy  performance 
prediction  programs  and  analyzes  the  characteristics  of  each.  This  information  should 
prove  a  valuable  contribution  to  future  fleet  Improvements  in  on-scene  prediction. 


C.  H.  Bassett 
Captain,  USN 
Commanding  Officer 
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M :  1 1  '  t«  :i:i!  n  nl  !  :k'I  i  <  ■ :  i  1  pn  it;  mills  \|  >1  1  I  ASDA,  AZ<  )|  :in<  I  S I  ’AM  ure  comparable, 
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:iii'l  miii|in  lealure-  "I  each  iHounim.  No  attempt  li:is  been  made  to  evaluate  the-  accuracv 
ot  output  values. 


In  an  open  environment  where  there  is  little  or  no  information  about  the  target,  the 
simulation  programs  give  varying  results.  The  data  in  this  report  show  that  these  programs 
evaluate  probability  of  detection  (P<>D1  differently  for  the  same  set  of  patterns  and  input 
values.  By  analvzing  cumulative  P(  >0  curves  there  are  general  relationships  that  can  be  ob¬ 
served  between  the  models:  (1)  TASDA  tends  to  be  conservative  in  comparison  with  AZ(  H 
and  SI’AIU  at  low  figures  of  merit  since  it  assigns  a  POD  of  0  or  1  during  a  Monte  Carlo  trial, 
(if)  SPAM  is  the  most  optimistic  of  the  three  at  low  figures  of  merit  since  it  uses  a  cumula¬ 
tive  time  window  and  a  weighted  average  to  obtain  POD.  (3)  At  all  frequencies  and  I’OMs, 
AZOI  shows  small  variance  between  probability  values  for  each  of  the  patterns  since  the 
calculation  of  POD  is  time  independent.  (4)  SPAM  is  more  conservative  in  its  POD  values 
at  high  I'OMs  since  the  cumulative  time  window  builds  detection  probabilities.  (5)  Both 
AZol  and  TASDA  have  a  sharp  broakolf  point  for  FOM  below  which  POD  values  decline 
rapiiliv  to  zero  because  ot  the  probability  of  detection  vs.  range  curves  used  in  each 
model.  Output  for  ADKPS  was  not  analyzed  since  the  number  of  patterns  it  evaluates  is 
limited  and  the  output  is  a  conditional  probability  of  detection.  Differences  between  com¬ 
puter  models  need  to  be  investigated  further.  However,  further  investigation  should  be 
based  on  a  real  world  application  of  the  sonobuoy  programs. 


I. 


LNTKODl'CTION 


Deciding  on  a  "best"  pattern  is  not  an  easy  task.  Ideally,  after  input  of  all  available 
data  and  execution  of  a  sonobuoy  program,  program  output  provides  an  optimal  sonobuoy 
pattern.  There  are,  however,  mnnv  complex  decisions  to  make  in  specifying  a  pattern 
including: 


where  to  drop  each  sonobuoy, 

number  of  sonobuovs  to  drop  in  a  pattern, 

sonobunv  depth  settings. 


length  of  time  to  monitor  a  pattern,  and 
number  of  patterns  to  drop. 


All  of  these  decisions  are  interrelated.  For  example,  the  number  of  patterns  depends  on 
the  size  of  the  search  area,  (lie  aircraft  on-station  time  and  the  availability  of  other  air¬ 
craft.  Additionally,  the  number  of  buoys  required  for  a  particular  pattern  depends  on 
sonobuoy  depth  settings  and  monitoring  time  for  each  buoy.  Operational  constraints  and 
mission  goals  must  also  be  considered  to  provide  the  decision  maker  with  the  best  possible 
choices  for  each  decision. 

Four  of  tiie  tactical  sonobuoy  computer  programs  which  aid  in  these  decisions  were  in¬ 
vestigated  to  determine  their  functional  effectiveness  as  an  aid  to  planning  search  patterns. 
The  programs  studied  were: 

1.  Tactical  Sonobuoy  Decision  Aid  (TASDA) 

2.  Algorithm  for  Zone  Optimization  and  Investigation  (AZOI) 

,1.  Search  Pattern  Assessment  Model  (SPAM) 

4.  Automated  Deployment  of  Sonobuoy s  (ADEPS) 

This  report  provides  a  description  of  the  simulation  techniques  used  in  each  model  in  Sec¬ 
tion  II,  Input  parameters  of  target  intelligence,  environmental  data,  buoy  data,  and  aircraft 
characteristics  are  compared  in  Section  III  followed  by  a  discussion  of  computer  output  from 
four  scenarios  in  Section  IV.  Sonobuoy  pattern  definition  is  included  in  Appendix  A.  The 
propagation  loss  environment  is  described  in  Appendix  B.  Output  data  analyzed  for  this  re¬ 
port  are  included  in  Appendix  C. 
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II. 


DESCRIPTION  OF  COMPUTER  SIMULATION  MODELS 


The  tactical  programs  ADEPS,  TASDA,  AZOI  and  SPAM  provide  measures  of  effective¬ 
ness  for  sonobuoy  patterns.  ADEPS  uses  a  mathematical  technique  based  on  the  symmetric 
properties  of  certain  sonobuoy  fields  to  calculate  conditional  probability  of  detection.  The 
other  models,  TASDA,  AZOI  and  SPAM,  use  Monte  Carlo  simulation  techniques.  In  these 
latter  three  models  satistical  distributions  are  used  to  simulate  initial  target  location.  If 
target  position  and  speed  can  be  estimated,  a  normal  statistical  distribution  is  used,  other¬ 
wise,  a  uniform  statistical  distribution  is  used  and  the  target  is  considered  equally  likely  to 
be  anywhere  in  a  search  area. 

Model  Description  for  ADEPS 

Two  search  methods  are  available  in  ADEPS  (Allison,  1970).  An  area  search  rou¬ 
tine  predicts  effectiveness  from  a  latticed  sonobuoy  pattern  with  equally  spaced  sonobuoys. 
This  routine  includes  the  following  possibilities  for  buoy  placement :  4-4,  2-2-2-2,  3-3-3, 
8-8,  5-6-5,  4-4-4-4,  6-6-6-G,  16-10,  11-10-11,  8-S-8-8.  A  barrier  search  routine 
assumes  that  a  submarine  is  on  a  course  perpendicular  to  a  row  of  sonobuoys.  Only  the 
area  search  method  is  examined  in  this  report. 

hi  the  ADEPS  area  search  routine,  symmetric  characteristics  of  the  sonobuoy  field 
are  used  to  divide  the  field  into  identical  subareas  as  shown  in  Figure  1.  A  lateral  range 
routine  computes  POD  versus  range  values  to  give  the  probability  that  a  sonobuoy  will 
detect  a  target  when  the  target  is  at  a  given  range  from  the  sonobuoy.  Individual  buoy 
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I  igure  1.  A  Symmetrical  Subsection  of  a  4-4- 1-4  Buoy  Field  in  ADEPS. 

Probability  of  Detection  is  calculated  lor  a  symmetrical  sub- 
area  of  the  buoy  field  and  is  used  as  indicator  of  effectiveness 
for  tile  entire  field. 

lateral  range  curves  are  combined  to  determine  probability  of  detection  at  equally  spaced 
points  ii.  a  subarea.  Conditional  probability  (P(t))  is  computed  from  the  sample  points  in  a 
.-.ubarea  where  Uie  combined  probability  of  detection  is  greater  than  or  equal  to  0.5.  The 
average  distance  D  between  consecutive  detections  is  computed  by  sampling  each  point 
sequentially  (as  shown  in  Figure  5t  and  averaging  the  distances  between  pnj)  >  0.5  sample 
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In  equation  1,  P(t)  is  a  conditional  probability  of  detection  in  a  subarea  during  time  t. 

P(t)  -l-e"T,t  Ttt  0  (  1) 

Although  it  is  not  the  probability  that  a  submarine  will  be  detected  bv  the  entire  sonobuov 
field,  it  is  a  measure  of  the  effectiveness  of  the  field.  The  final  product  of  the  model  is 
the  average  conditional  probability  of  detection  (P)  given  that  a  target  enters  a  symmetrical 
subsection  of  the  buoy  field,  ft  is  determined  by  averaging  the  values  for  P(t)  over  the  on 
station  time  period. 

Model  Description  for  TASDA 

TASDA  operations  (NADC,  197(1)  are  divided  into  two  functional  objectives:  (1)  crea¬ 
tion  of  a  sonobuov  tactics  file,  and  (2)  application  of  Monte  Carlo  simulation  to  produce 
optimum  sonobuov  field  deployment  for  ASW  threats.  The  tactics  file  is  created  by  a 
geometry  definition  program  (GEOMT).  TASDA  uses  Monte  Carlo  game  theory  to  stimu¬ 
late  100  random  tracks  of  a  target  through  a  sonobuoy  field  selected  by  GEOMT.  Aircraft 
time  on  station  is  divided  into  time  steps.  The  position  along  a  target  track  is  computed 
for  each  time  step.  If  signal  excess  is  greater  than  zero  for  two  consecutive  time  steps, 
a  detection  is  recorded.  The  frequency  of  detected  target  tracks  to  total  tracks  is  used 
as  an  estimate  of  the  expected  and  cumulative  probability  of  detection.  The  total  time 
a  target  is  detected  bv  one  or  more  sonobuoys  divided  by  the  number  of  unbroken  detec¬ 
tion  intervals  is  output  as  mean  holding  time.  The  time  of  first  detection  for  each  target 
track  is  recorded  and  the  total  of  these  times  divided  by  the  number  of  detected  tracks 
is  output  as  mean  time  to  first  detection. 


Model  Description  for  A/.<  d 


A/.i  >1  (Birnbaum  and  Dorav,  1!)7'>1  is  ;i  taclics  package  which  includes  four  subprograms 
1  .<  ><  >P,  oPTI,  WAI  T,  and  CAM!  .  Ixpected  probabi lit v  of  detection  <  I- POD)  bv  a  sonobuoy 
field  is  computed  in  I.OOl*.  OPTI  determines  the  pattern  spacing  which  provides  maximum 
I  Poll.  W Al  l  calculations  include  cumulative  probability  ol  detection.  CAM!  uses  Kl*o|) 
obtained  1>\  !  <  >i  )|>  to  estimate  an  optimum  buoy  spacing  for  a  given  pattern  with  several 
different  sonobuoy  chnraet eristics,  such  as  figure  of  merit  (FOMi  or  buoy  depth. 


In  A'/.nl,  I  pol)  is  estimated  by  Monte  Carlo  simulation.  A  lateral  ratine  curve  of 
probability  of  detection  versus  signal  excess  is  given  bv  a  receiver  operating  curve  (TfOC't 
which  is  determined  irom  a  probability  ol  false  alarm  ol  1<C^  and  a  specific  spectrum 
analv/.er  set  lor  a  live-minute  integration  lime.  Signal  excess  versus  range  is  deter¬ 
mined  from  the  F<  *M  and  the  propagation  loss  curve.  The  two  tables  of  values  are  combined 
to  produce  Pol)  v  ersus  range  values.  H.ange  from  a  sonobuoy  is  determined  from  a  statis¬ 
tical  target  distribution.  Detection  probability,  Pi,]*  for  the  jt h  sonobuoy  at  range  rj  from 
the  target  is  found  from  POD  versus  range  values  (Figure  .Mi.  In  equation  2,  Z(Xj,  Yj)  is  the 
expected  value  for  detection  when  the  target  is  at  a  partieidar  position  (Xj,  Yj).  The  sub¬ 
script  i  is  an  index  used  to  number  target  positions  from  1  to  N. 
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n  =  number  of  sonobuoys 


KPOD  for  a  pattern  is  calculated  by  talcing  a  random  sample  ol  N  target  positions.  Proba¬ 
bility  is  calculated  for  each  position  i  as  in  equation  2.  The  probabilities  are  summed  and 
divided  by  X  to  calculate  KPOD  in  equation  M. 

lPOf  =  z( x  ,y  )  *  ?.•  >:  ,v  '  <■  j- x  :<  }  + . z{>:..,y.1)  /n 

H  i  ' 


N  =  number  of  target  positions 

Cumulative  probability  of  detection  (CPODi  i-  also  estimatetl  by  Monte  Carlo  simu¬ 
lation.  The  mathematical  function  used  in  live  computation  ol'  CPOD  is  binomial  and  gives 
tile  probability  oi  at  least  one  detection  in  X  trials.  A  time  tj  is  chosen  randomly  from  a 
time  interval  I  using  a  random  number  generator.  The  interval  I  is  defined  by  equation  A 


T  -  dir  craft  time  iute 
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Ar  i.ini  er'ainty  about  the  target  posi'io.i  at  time  tj  is  caleul  tied  from  a  normal  statistical 
distribution.  Detection  probability  at  a  position  i.-  ■.  aiculn.ed  as  in  equation  2.  The  maxi¬ 
mum  expected  value  fm  detection  bv  sonobuoy  i  (Figure  D  over  X  target  tracks  at  each 
time  tj  is  summed  to  yield  cumulative  probability  of  detection  in  equation 
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The  CPOD  in  equation  5  is  an  estimate  for  the  lower  limit  of  cumulative  probability  of  de¬ 
tection  at  time  t.  A  large  number  of  target  tracks  must  be  generated  to  obtain  an  estimate 
of  CPOD  with  a  small  variance.  The  number  of  simulations  required  depends  on  the  signal 
excess  and  the  desired  accuracy  of  output  data.  The  probability  of  detecting  a  target  that 
has  not  yet  been  detected  is  calculated  from  CPOD  using  Bayes  theorem  of  probability. 


Figure  3.  A  Graph  of  POD  vs.  Bange  Values  Computed  in  the  AZOI  Tactics  Package. 


Figure  4.  Maximum  POD  Pjj*  of  3  Target  Tracks  at  Time  tj  for  a  Sonobuoy  j 
in  AZOI. 
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Model  Description  for  SPAM 


The  SPAM  program  computes  (''POP  bv  Monte  Carlo  simulation  (Marin,  1979).  The 
range  of  the  target  from  a  sonobuov  is  found  bv  a  stochastic  process.  Signal  excess  (SE) 
at  time  t  is  found  from  a  normal  distribution  with  mean  SK(r)  and  a  standard  deviation 
of  li,  >,  10,  or  12  tlB.  The  probability  that  the  signal  excess  is  greater  than  zero  at  range 
r  is  given  by  equation  (»  (Marin,  197(>). 


,  r  (x  - 

P  ISEtr.)  0 1  =  - ---  /  e  -  - ---  dx  (6 

r  ■  1  j  2 

-  -  standard  deviation 
of  signal  excess 

To  vise  equation  (>  in  a  simulation  model,  it  is  necessary  to  introduce  a  time  element.  The 
expression  Pj(tji  is  introduced  as  the  probability  that  a  target  is  not  detected  at  riuige  rj 
during  time  t;.  The  following  equations  are  assumed  for  the  time  element. 


iT\t.)  "  ..  t .  +  0  ( t .  )  where  lim  0  ( t }  / 1  -  0  (7a 

J  1  t  •  0 

P.(t.)  +  M  j  ( t  • }  -  1  where  (E (0)  =  1  (71 

^jl'tj  +  t)  •-  *  m . ( ■  1 1  ;  7c 

liquation  7a  states  that  the  probability  of  detection  Pj(tj)  is  proportional  to  some  constant 
A  -  plus  some  higher  order  terms.  The  higher  order  terms  in  equation  7a  can  be  neglected 
for  a  small  time  period  t.  Equation  7b  states  that  the  target  is  either  detected  or  not 
detected  where  the  probability  Mj(0)  of  no  detection  at  any  range  with  no  time  elapsed  at 
that  range  is  one.  Equation  7e  assumes  that  no  detection  at  time  t  and  no  detection  at  time 
t-  At  are  independent  events.  The  assumptions  lead  to  an  exponential  statistical  distribu¬ 
tion  for  Pj(tj>.  The  SPAM  program  uses  the  equation  (Marin,  197-~>) 

P.it. .  )  -  1  -  e  for  t.  0  *g9 

J 

in  a  simulation  technique  called  the  cumulative  time  window  (CTW).  The  length  of  time  t() 
necessary  for  a  sonobuov  j  to  detect  a  target  and  the  length  time  t^  between  time  steps 
are  input.  Probabilities  Pjdjl  are  computed  and  compared  with  random  numbers  to 
determine  when  a  detection  occurs.  A  program  trial  ends  either  when  a  detection  occurs 
or  when  the  aircraft  on  station  time  expires.  The  upper  limit  for  cumulative  probability  of 
detection  is  obtained  In  combining  single  buoy  probabilities  independently  in  equation  9. 


rP:-,r',„, Af , ;  ■  -  i  -  Kfj)  .  .  i:  -  ...c-  puy;  (9) 

n  -  number  of  sonobuoys 

A  lowci  limit  for  probability  oi  detection  at  time  tj  is  obtained  assuming  complete 
dependence  as  in  equation  PE 


i1AX(f-.(ti)) 


(  POD  . 


mi  n 


1, 


n  =  number  of  sonobuoys 

SPAM  outputs  cumulative  probability  slightly  less  thiui  hallway  between  the  two  extremes 
using  equation  11. 


(IPOD  =  .45  ( CPOD  +  CPOD  ,  ) 
max  min 


CPOr 


mi  n 


Ml 


Summary  of  Model  Descriptions 

Knell  of  the  models  calculates  probability  of  detection  at  some  time  tj  in  order  to 
arrive  at  an  estimate  for  KPOD.  ADKPS  calculates  the  average  conditional  probability 
P  in  a  symmetrical  subarea  of  a  buoy  pattern.  TASDA  computes  probability  of  detec¬ 
tion  bv  one  or  more  buoys  (Pl)lt  using  Monte  Carlo  simulation  of  100  target  tracks  anti  a 
detect  i  POD  -  D  or  no  detect  (POD  -  Oi  analysis  for  ea^h  track.  A'/.oi  uses  a  binomial 
distribution  function  to  calculate  KPoD  where  a  value  along  a  time  track  is  chosen  random- 
lv.  SPAM  uses  Monte  Carlo  simulation  and  averages  Pj(tj)  (previously  defined")  obtained  by 
assuming  dependence  of  sonolnioys  with  Pj ftj )  obtained  by  assuming  independence  of  sono- 
buovs. 

Inch  of  the  models  estimates  cumulative  probability  of  detection.  ADKPS  uses  an 
extvinential  expression  to  estimate  CPOD.  TASDA  uses  a  detect -no  detect  Monte  Carlo 
mathematical  technique.  AZoi  uses  some  of  the  KPOD  calculations  to  obtain  an  estimate 
for  the  lower  limit  of  CPOD.  In  the  uniform  case  CPOD  is  time  independent  and  is  equal 
to  KPoD.  SPAM  uses  a  cumulative  time  window  (CTW)  defined  by  tQ  and  to  calculate 
CPOD. 


Ill. 


SCI- NAT? IO  DKVINITH  >N 


Five  sonobuoy  patterns  employed  by  ASW  tacticians  were  chosen  lor  data  ;mal\sis  in 
scenario  1:  the  circle  (ellipse  in  the  normal  ease),  1-1  1-1,  chevron  (chevron  skew  in 
the  normal  ease),  3-0-3,  and  brushlae  patterns.  Their  coni' Lgu rations  are  dclincd  in 
Appendix  A.  The  same  input  lor  target  intelligence,  environmental  data,  and  aircraft 
characteristics  were  used  so  as  to  evaluate  the  effectiveness  of  each  pattern.  Output 
was  generated  for  one-half  hour,  one  hour,  and  four  hours  of  flight  time  in  all  cases 
preceding  the  aircraft's  arrival  on-station.  The  aircraft  remained  on-station  four  hours 
in  all  cases  for  scenario  1.  The  circle,  ellipse,  1-4— 1-4  and  chevron  skew  patterns 
were  evaluated  in  scenario  2.  The  on-station  time  was  extended  from  lour  to  l(i  hours 
in  the  second  scenario  and  output  was  generated  for  one  hour  time  late,  hi  scenario  .4, 
the  circle,  ellipse,  4— 1-4-4  and  chevron  skew  patterns  were  searched  in  quarters  In 
four  planes  with  one  hour  time  late  and  tour  hours  on  station.  Scenario  4  used  the 
same  target  data,  buoy  data,  and  aircraft  data  input  as  scenario  3,  however,  environ¬ 
mental  input  was  changed  for  oilier  comparison  purposes. 

Scenarios  1-3  were  based  on  the  following  environmental  input.  An  XBT  temperature 
profile  to  730  meters  recorded  0  July  1971  in  tlie  Pacific  Ocean  at  latitude  33°  3s'  North 
and  longitude  123°  39'  West  was  input  into  the  Integrated  Command  ASW  Prediction  System 
(ICAPS)  PROFGKN  computer  program.  This  program  merged  the  observed  depth- 
temperature  pairs  with  a  deep  historical  profile  to  yield  a  total  temperature  profile  with 
bottom  depth  at  4000  motors.  The  program  then  calculated  a  sound  speed  prolile  from  the 
total  temperature  profile  and  historical  salinity  values.  The  sound  speed  profile,  which 
exhibited  a  depth  excess  indicative  of  convergence  zone  environment  (Figure  3),  was  input 
into  the  ICAPS  Fast  Asymptotic  Coherent  Transmission  (FACT)  program  to  determine 
acoustic  propagation  loss.  By  varying  die  target  depths,  receiver  depths,  anti  detection 
frequencies  inserted  into  the  FACT  program  nine  propagation  loss  (PL)  curves  were 
generated.  These  curves  appear  in  Appendix  B. 

The  environmental  input  for  scenario  I  was  a  Gulf  Stream  profile  from  latitude  39°  3d' 
North  longitude  70°  29'  West  recorded  17  November  1909.  If  was  merged  with  ICAPS  deep) 
history  data  to  produce  a  total  profile  with  bottom  depth  2400  meters  (Figure  (5).  An  At¬ 
lantic  slope  water  profile  at  latitude  39°  33'  North,  longitude  71°  30'  West  recorded  18  Nov¬ 
ember  1909  was  merged  to  produce  a  total  profile  with  the  same  bottom  depth  (Figure  7i. 
Two  different  profiles  were  used  in  scenario  1  in  order  to  show  the  effect  of  varving  the 
environment  on  the  tactical  programs.  Since  scenario  4  uses  environmental  input  which 
is  different  from  the  other  •'!  scenarios,  its  output  does  not  directly  compare  with  output 
from  the  first  three  scenarios. 
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Figure  f>.  Atlantic  Culf  stream  sound  speed  profile  used  to  generate  propagation 
loss  curves  for  scenario  4.  (The  predominant  mode  of  sound  propagn 
tion  i  s  bv  tlircct  path. ) 


UELOCITY  SOUND  UELOCITY  PROFILE  LAT  3936N  L 

UELOCITY  n/SEC  DATE  171169 
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All  tour  tactical  sonobuoy  computer  programs  require  buoy  input,  target  input,  and 
environmental  input.  ADEPS  does  not  require  aircraft  input.  Some  input  values  are  preset 
in  each  program  so  that  the  user  need  not  specify  a  long  list  of  values.  Data  processing 
inputs  for  the  four  scenarios  considered  in  section  V  are  given  in  Table  1. 

Buoy  Input 

Buoy  input  which  define  sonobuov  patterns  for  evaluation  is  shown  in  Table  2.  In 
ADEPS  the  total  number  of  sonobuoy s  input  must  be  a  multiple  of  eight.  Area  dimensions 
are  required  and  are  used  to  calculate  a  buoy  spacing.  The  buoys  are  equally  spaced  to 
cover  an  entire  area.  ADEPS  is  die  most  limited  program  in  terms  of  buoy  input,  but  the 
symmetric  patterns  allow  rapid  processing.  The  only  geometries  available  in  die  ADEPS 
preset  pattern  package,  applicable  for  Uiis  study,  were  the  4-4-4— 4  pattern  and  the  5-6-"> 
pattern  with  a  uniform  target  distribution. 

TASDA  requires  the  number  of  sonobuoy s,  buoy  position  assignment,  the  number  of 
different  spacings,  and  buoy  spacing  limits  as  input.  Patterns  are  input  using  a  sonobuoy 
position  planner  chart.  Optionally,  die  user  may  input  buoy  positions  in  x-y  coordinates  to 
define  a  pattern.  GEOMT  is  a  preprocessor  program  for  TASDA  and  creates  a  tactics  file 
of  up  to  20  sonobuoy  patterns  with  as  many  as  64  buoys  in  each  pattern.  This  preprocessing 
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N  .'iHitc  combines  sio  cases  shown  in  Table  1  info  11>2  runs  for  processing.  Hun  time  and 
I'tiov  inpu!  time  are  reduced  since  buov  positions  are  simple  recalled  from  OEOMT  rather 
t  hail  rout ; nel\  entered. 


Table  2.  Buoy  Input 


ADEPS 

TASDA 

AZOI 

SPAM 

Coordinates 

X 

X 

X 

Spacing 

X 

X 

X 

X 

Buoy  fi  ansl'ormation 

X 

X 

Preset  Patterns 

X 

Geometries 

X 

Ambient  Noise  SDev 

X 

Di  AZ OI  a  pattern  geometry  is  defined  by  specifying  a  set  of  x-y  coordinates  for  each 
buoy.  A  lactor  for  converting  relative  distance  to  nautical  miles  must  be  specified.  Four 
ty  pes  of  transformations  may  change  the  basic  pattern:  x-spacing  expansion,  y-spacing 
expansion,  pattern  shift  along  the  x-axis,  and  rotation  of  the  pattern  about  the  origin. 

In  SPAM  the  user  specifies  the  number  of  sonobuoys  mid  the  x-y  coordinates  in  nauti¬ 
cal  miles  for  each  buoy.  Buoy  coordinates  are  input  in  the  order  in  which  the  buoys  are 
deployed.  As  the  buoys  are  dropped,  the  user- specified  radio-frequency  (RF)  range  is 
used  to  determine  if  a  sonobuoy  can  be  monitored.  After  the  pattern  is  complete,  SPAM 
simulates  monitoring  constraints.  The  user  specifies  the  average  number  of  buoys  that 
an  aircraft  can  monitor  and  the  monitoring  time.  Using  the  value  for  the  average  number 
of  buoys,  the  program  randomly  selects  buoys  from  the  pattern  and  only  these  buoys  will 
be  used  to  determine  if  a  detection  occurs. 

Target  Input 

There  are  twenty  target  inputs  as  shown  in  Table  3.  In  ADEPS  only  four  of  the  inputs 
are  applicable.  Target  speed  is  input.  The  area  search  option  corresponds  to  an  evenlv 
spaced  position  input  where  the  target  remains  hi  the  search  area  at  fixed  positions.  The 
barrier  search  option  corresponds  to  the  univariate  normal  input  where  the  target  is 
assumed  to  progress  towards  a  barrier  with  a  target  track  perpendicular  to  the  barrier. 
Only  the  evenly  spaced  input  is  used  in  this  study  to  search  an  area  with  and  1-1- 1-1 

patterns . 

Twelve  of  the  twenty  target  inputs  apply  to  TASDA.  Target  speed  is  required.  The 
target  type  must  be  specified  as  nuclear  or  conventional,  holding  or  transiting.  A  nuclear 
holding  target  was  used  for  this  study  (see  Table  1).  Target  location  may  lie  described 
in  two  ways: 

1.  A  uniform  time-dependent  density  function  assumes  that  the  target  density  is 
random  over  a  designated  rectangle  (Figure  X). 
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2.  A  Bessel  normal  distribution  function  is  generated  from  an  initial  bivariate 
normal  density  function  that  is  expanded  at  a  rate  equal  to  input  target  speed 
(Figure  !)). 

Both  target  location  functions  were  used  in  this  study. 


Table  3.  Target  Input 


Target  Movement 

ADEPS 

TASDA 

— 1  - 

AZOI 

SPAM 

Course 

X 

X 

Course  Standard  Deviation 

X 

Course  Limits 

X 

Speed 

X 

X 

X 

X 

Speed  Standard  Deviation 

X 

X 

Velocity  Vector 

X 

Conventional 

X 

Nuclear 

X 

Holding 

X 

I  r  tut  si  ting 

X 

Distance  Traveled 

X 

Snorkel  Cycle 

X 

Target  l  ocation  standard  Deviation 

X 

X 

Initial  Target  Location  Density  Functions 

Kvcnlv  spaced  positions 

X 

t  ni  form  Time -dependent 

X 

l  nifoini  Time-independent 

X 

X 

Tain  ariatc  normal 

X 

X  ! 

!>•*.'  ^ci  I’.unruJ 

X 

X 

Hi',  ariatc  normal 

X 

X 

Time-dependent  bivariate  normal 

X 

i 

J 

\X<  d  uses  seven  ol  lwent\  target  inputs.  T;u'get  speed  and  speed  standard  deviation 
are  input  di  recti \  or  m  term.-  ot  ,i  speed  covariance  matrix.  I  he  following  four  target  lot  a- 
lion  ’ 1  •  :  a  functions  :ot  used  to  predict  initial  target  location  for  Monte  C  arlo  simulation: 


RANDOM  TARGET  TRACK 


Figure  s.  Fniformly  Distributed  Initial  Target  Area  in  TASDA.  The 

target  location  is  random  over  a  designated  area.  The  same 
type  of  target  area  definition  is  available  in  AZOI  and  SPAM. 
However,  the  statistical  distribution  is  time-dependent  for 
TASDA. 


Figure  9.  Normally  Distributed  Initial  Target  Area  in  TASDA.  The  ini¬ 
tial  target  location  is  normally  distributed  around  (XPT,  YPT) 
with  standard  deviation  SGMAX  in  the  x-direction  and  SGMAY 
in  the  y  direction.  The  same  type  of  area  definition  is  avail¬ 
able  in  AZOI  and  SPAM.  However,  the  statistical  distribution 
is  bivariate  normal  for  SPAM. 


1.  A  uniform  time-indei>endont  density  function  assumes  that  target  density  is  con- 
st;uit  over  a  designated  area. 

2.  A  Bessel  normal  density  function  is  generated  from  an  initial  bivariate  normal 
density  function  and  expands  at  a  rate  equal  to  the  best  estimate  of  target  velocity 

:>.  A  time -dependent  bivariate  normal  density  function  predicts  a  target  position 

when  target  fracking  has  been  successful.  This  function  describes  a  target  being 
tracked  when  it  is  imperative  to  re-establish  target  contact. 

i.  A  bivariate  normal  distribution  function  assumes  that  a  target  is  normally  dis¬ 
tributed  in  an  area  generated  from  a  single  DII'AH  buoy  fix,  a  SOSl’S  contact,  or 
radar  contact.  This  distribution  does  not  expand  with  time. 

The  first  and  second  functions  were  used  for  uniform  and  normal  cases  in  this  study 
(Table  1). 

SPAM  uses  seven  of  the  twenty  inputs  listed  in  Table  3.  Required  input  is  speed, 
speed  standard  deviation,  course  and  course  standard  deviation.  Target  location  can  be 
assigned  using  three  distribution  functions: 

1.  A  uniform  time -independent  distribution  function  with  the  property  that  a  target 
is  equally  likely  to  be  at  any  point  within  a  bearing  box  of  length  L.  and  half¬ 
width  H . 

2.  A  bivariate  normal  distribution  function  which  has  a  user  specified  mean  and 
standard  deviation. 

;i.  A  univariate  norm  til  distribution  function  which  has  a  normal  distribution  in  an 
\  direction  and  a  uniform  distribution  in  a  y  direction. 

function.-,  1  and  2  were  used  for  data  analysis  in  this  study. 

Ae<  uslic  Input 

Acoustic  input  for  each  program  is  listed  in  Table  -J.  The  acoustic  input,  source 
level  (SI.),  ambient  noise  (AX),  and  recognition  differential  (RD)  are  used  to  determine 
R0T  .is  in  equation  12. 

run  !.l  -  -  rl-  (12; 


Probability  of  detection  is  defined  to  be  < ) .  A  when  signal  excess  (SE)  in  equation  13  is  zero 
for  propag  ition  loss  (PI.)  at  range  R. 

!  =  fof  -  n  i :: 


p. 


Ln  each  program  a  probability  curve  is  used  to  assign  probability  values  to  values  of  signal 
excess.  When  signal  excess  is  greater  than  zero,  the  probability  values  will  be  greater 
Uian  n..j.  When  signal  excess  is  less  than  zero  probability  values  will  be  less  than  0.5. 


Table  L  Acoustic  Input 


AD  EPS 

TASDA 

AZOI 

SPAM 

ITop  l  oss  (PL) 

X 

X 

X 

X 

S'  'ii  re  1  i'v<  1 

X 

Mubieni  Noise 

X 

Kecognition  Ditierentkil 

X 

Figure  of  Merit 

X 

X 

X 

Multiple  FO.M 

X 

X 

X 

K<  >C  Curve 

X 

Multiple  PL 

X 

X 

X 

Convergence  Zones 

X 

X 

ADI  PS  considers  propagation  loss  (PI.)  curves,  source  level  (SL),  ambient  noise  (AN), 
and  recognition  differential  (HD),  but  ADKI’S  does  not  have  a  provision  for  time  late.  As 
a  result  ADI  PS  could  only  be  applied  for  54  of  the  810  possible  cases.  There  is  no  pro¬ 
vision  for  multiple  frequencies  or  multiple  FOMs  and  a  separate  run  was  required  for  each 
case.  In  order  to  reduce  the  number  of  runs,  an  FOM  of  75  was  selected  with  a  frequency 
of  500  Ilz,  Pighteen  computer  runs  were  processed  for  scenario  1. 

The  TASDA  program  cun  accept  up  to  three  propagation  loss  curves  and  up  to  four 
FOMs.  Those  two  features  reduced  the  number  of  runs  required  from  102  to  18.  The  user 
may  input  convergence  zone  ranges  and  widths  instead  of  propagation  loss  curves. 

In  AZOF,  a  receiver  operating  characteristic  (HOC)  curve  is  required  input.  Ibis 
curve  gives  probability  values  as  a  function  of  signal  excess.  Other  inputs  are  a  table  of 
propagation  loss  versus  range  and  a  figure  of  merit.  Two  subprograms  of  AZ()1  were 
modified  in  Ma\  1!)7H  to  process  multiple  FOMs,  multiple  propagation  loss  curves,  and 
multiple  patterns.  The  number  of  required  data  sets  for  input  were  therein-  reduced  from 
s  10  to  eight . 


-  i 1  \  \  |  accept  s  i  ii u-  propac.ii  inn  lit-.-  Ii !  ■  a  1 1 !'  1 1  <  • ! . i < ■  i ic \  and  i  eeei\  e  r-  sou  rc  i  1 1 1  (  li  - 
si;t  .  cit  pih  idem  i  l  ien!  pm  immi»  r-.  Am  immln  r  nl  pr<  >  p :  i  lj.  :  i  t  ion  loss  curve-  in.r  be  input 
to  tin  lib- william  nimilu  i'  oi  li'upi  m  ii  !  '  •  'lo.inim  uutoniat  ienll  \  -enrelie-  t  Ii  rough 
i Ik  input  propagation  loss  till-  ...  v  ill..  a|»propriat «  pr< vacation  loss  aini.  1  !w  main  pro¬ 
gram  SI’AM  u  is  modified  in  .srpltmla-r  luTT  for  multiple-  l  <>M,  multiple-  time-  late,  and 
I’v.iltiple  propagation  loss  input.  I  he  numl't-r  of  executions  was  reduced  from  "10  to  ei o lit . 

\  pi  opap.U ion  loss  cum  is  cins-ilied  as  -hallo" . --hnl lov  (S-Sl  lor  a  shallow  target  and 
r,  cciu  i',  il< -  p-det  p  i  |)-|)i  lot  a  deep  ’  •  reel  and  receiver,  or  shallow-deep  t S—  I >)  for  either 
a  -hallo. target  or  receiver  or,  output. 

Airt  raft  Input 

A  DHL'S  docs  not  consider  aircraft  parameters  except  tor  monitor  degradation  as  shown 
in  Table  a.  For  TASDA,  tile  flight  path  is  described  onl\  to  Die  extent  of  specif.!  ing  die 
■  •lapsed  lime  before  arrival  on  station.  The  radio  Iroquonci  (RF)  nuigc  is  the  maximum 
allowable  distance  between  Uic  aircraft  and  a  sonobuoy  which  permits  radio  reception, 
lor  TASDA.  the  position  of  the  aircraft  is  fixed  at  the  center  of  the  field  for  tile  simulation. 
\/*  'I  uses  time  late,  aircraft  velocity,  and  on-station  time  to  calculate  I'l’nD.  SI’AM  in¬ 
cludes  time  late  and  time  on  station  as  well  as  aircraft  "round  speed.  The  user  specified 
around  -peed  in  SI’AM  determines  the  time  at  which  each  buoy  comes  up.  Onl\  those  lnio\  s 
.  i  thin  a  a  set- specified  l!F  ranee  are  monitored.  After  all  buoys  are  in  the  water,  die  user 
cai!  specif'  -i  sonobuoi  monitoring  scheme  that  reflects  RF  constraints. 


Table  5. 

Aircraft  input 

1 

-  ADI.  1>S 

TASDA 

A  7.01  [  SI’AM 

1  i ;  I  1  Tinge 

X 

1  -X 

1  lime  Late 

X 

1 

i 

X  .X 

1  Tin. ■  -  on  Station 

X  ,  X 

'.lu-iitornig  Sequence 

j 

1  X 

•  i  -i.iv!  Spiced 

1  x 

Monitor  Degradation 

1 

1 

!  x 

_ L 

_J _ 

1 

— — — 1  —  I.  — 

Sui i a  Mo; tel  Input 

i.a-  ■  oi  '.la  pi  op  nun.-  uses  dillcrent  mathematical  loi  nuilatioiis  to  describe  a  tactical 
-it.  alio  c.  ii  ■•pen  ocean  einironmci'l.  Die  lormulativum.  are  based  on  assumptions  about 
tile  \  annum  ■  -t  die  J  Cu.M.  Ui«  amount  oi  signal  e.\Cess  necessary  to  detect  a  target,  ;uul  die 
-‘cpendi--  c  e  (or  independence)  ol  sonobuoy.-  collectivc-iv  contributing  to  a  detection. 


'.!)!  1  i'U'M'i  puff urns.  Inc  preset  package  cannot  be  modified.  TASDA  also 

"ive^ct  pat le rn  package,  but  t he  package  can  be  mollified.  AZoi  and  SI’AM  have  no 
preset  pattern  files.  Al)l  PS  has  a  normal  initial  target  location  distribution.  TASDA, 

Y/i  >!  and  SPAM  all  have  options  for  uniform  target  location  distributions  and  normal  loca- 
■  inn  dist  ri  but  ions. 

\D  models  input  a  propagation  loss  curve.  TASDA  may  input  propagation  loss  curves 
at  Po-ee  different  depths.  Therefore,  one  sonobuoy  pattern  may  have  sonobuoys  placed  at 
i!i tb.- rent  depths.  AZ(  '1  inputs  a  receiver  operating  characteristic  (HOC)  curve  which  assigns 
■"  t-robability  to  v  alues  of  SI  .  In  ADKPS  and  SPAM  this  curve  has  a  normal  density  distri¬ 
ct  ion  with  probability  e<p: rd  to  0.5  at  the  mean  of  the  normal  distribution.  In  TASDA 
cumulative  probability  i,->  determined  from  a  Pieian  density  Distribution.  TASDA  assumes 
till  ^.uioouoys  lire  activated  at  the  same  time.  In  SPAM,  the  user  specified  ground  speed 
detei  mines  the  time  each  buoy  comes  up. 

V.  ANALYSIS  ( >r  HKSI't.TS 

Two  types  of  program  out  pi'1  are  compared  by  this  study:  expected  probability  of  de- 
:  ection  i  PP<  )Dt  and  cumulative  probability  of  detection  (CPOD).  F.POD  follows  the  statistical 
definitions  for  expectation.  CIYiD  is  computed  from  a  probability  density  function  which 
follows  the  statistical  definition  of  a  density  function.  Each  program  uses  different  assump¬ 
tions  to  arrive  at  EPOD  and  different  density  functions  to  calculate  CPOD.  Since  methods 
of  calculations  are  different  in  each  program  the  output  have  different  numerical  values. 

The  analysis  in  this  section  examines  mathematical  techniques  to  explain  the  CPOD  and 
EP()D  numerical  output. 

A  list  of  all  program  output  is  given  in  Table  6.  Although  CPOD  and  EPOD  para¬ 
meters  in  Table  (i  are  listed  for  TASDA,  AZOI,  and  SPAM,  numerical  output  for  these 
measures  of  pattern  effectiveness  are  not  the  same.  The  output  values  were  analyzed  in 
reference  to  four  scenarios.  Three  target-receiver  depth  combinations,  shallow-shallow 
tS— S'),  shallow-deep  (S-D),  and  deep-deep  (D-D),  were  used  to  generate  propagation  loss 
cinnxs  for  the  first  3  scenarios.  Probability  of  detection  was  predicted  for  three  FOMs  and 
for  a  low  frequency  range,  middle  frequency  range,  and  a  high  frequency  range  in  scenario 
1,  2  and  3.  Scenario  1  uses  two  target-receiver  depth  combinations,  two  frequencies  and 
one  I'OM. 

Scenario  1 

Expected  probability  of  detection  was  calculated  for  five  patterns  with  a  uniform  and 
normal  initial  target  location  distribution.  Time  late  was  not  a  factor  in  the  uniform  case 
because  that  parameter  did  not  change  the  initial  target  location  distribution  for  AZOI  and 
SPAM.  Therefore,  values  for  1  and  4  hrs  time  late  have  been  omitted  and  only  values  for 
30  min  time  late  appear  in  Appendix  C.  Cumulative  probability  of  detection  was  analyzed 
using  the  method  shown  in  Figure  C-l.  Csing  this  method  cumulative  probability  of  detec¬ 
tion  was  converted  to  a  single  eventual  probability  of  detection  value  for  comparison  pur¬ 
poses.  Both  the  eventual  POD  values  and  the  original  CPOD  curves  appear  in  Appendix  C 
for  TASDA,  AZOI,  and  SPAM.  The  18  CPOD  curves  that  were  generated  for  ADEPS  are 
shown  in  Figure  10  with  the  corresponding  curves  for  TASDA,  AZOI  and  SPAM.  Since 
ADEPS  could  process  only  two  of  the  patterns  the  results  were  not  analyzed  for  inclusion 
in  Table  7.  However,  the  relationship  between  ADEPS  output  and  output  from  the  other 

iy 


three  models  can  be  deduced  from  Figure  10.  one  characteristic  of  the  ADFPS  curves  is 
lhat  tliev  alwavs  rise  to  100'  I >  1 ) . 


Table  G.  Measures  of  Effectiveness 


ADK  PS 

TASDA 

AZOI 

SPAM 

PROGRAM  OUTPUT 

X 

X 

X 

CPOD 

X 

X 

X 

EPOD 

X 

PD2 

X 

PD3 

X 

MIIT1 

X 

MHT2 

X 

MHT3 

X 

X 

MHFD 

/\ 

p 

X 

P(t) 

X 

SIGMA 

_ 

* _ 1 

_ 

P(DE/NDT) 

C  POD 

I.  POD 

PD2 

PD3 

MIIT1 

MIIT2 

MIIT3 

MTFD 

A 


lJ(t> 


STOMA 

P(DL/XDT) 


Oniform  <  ase 


cumulative  probability’  as  a  function  of  time  on  station 
expected  probability  of  detection  on  one  or  more  sonobuovs 
probability  of  detection  on  two  or  more  sonobuoys 
probability  of  detection  on  three  or  more  sonobuoys 
mean  holding  time  on  one  or  more  sonobuoys 
moan  holding  time  on  two  or  more  sonobuoys 
merm  holding  time  on  three  or  more  sonobuoys 
mean  time  to  first  detection 

average  conditional  probability  in  a  section  where  die  submarine 
is  assumed  to  be  a  random  target 

cumulative  probability  that  a  submarine  has  entered  part  of  a 
buoy  field  where  signal  excess  SE  is  greater  than  zero 
standard  deviation  of  POD 

probability  dial  a  target  will  be  detected  given  that  it  has  not 
been  detected  in  die  past 


r;p<  d)  \ alues  for  the  uniform  case  were  averaged  over  three  FOMs  and  three  frequen¬ 
cies  and  the  averages  are  present  in  Table  7.  For  a  shallow  target  and  receiver,  TASDA 
diuws  P*.c  1-1  — i— I  pattern  as  most  effective  over  the  frequencies  and  FOMs  with  a  value 
if  l.o  |  I'op,  For  the  shallow-deep  target-receiver  combination,  TASDA  shows  the 
chevron  pattern  as  most  effective.  For  a  deep-deep  target  receiver  propagation  loss 
curve,  TASDA  selected  the  chevron  pattern,  A r/,(  n  selected  the  circle  pattern  and  SPAM 


ADI';  i 


•  IASDA 
AAZOI 
OS  I ’AM 

5-6-5  PATTERN 


Figure  10.  Probability  of  Detection  versus  Flapsecl  Time  on  Station  for 
Scenario  1  with  a  Uniform  T  arget  Location  Distribution. 


21 
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selected  the  4-4-4-4  pattern  as  most  effective.  Note  that  there  is  very  little  agreement  in 
pattern  selection.  This  is  because  of  the  different  modeling  in  each  program.  Since  the 
patterns  were  defined  at  a  near  optimum  sonobuoy  spacing,  the  models  give  about  the  same 
probability  averages  for  each  pattern  they  evaluate.  However,  the  difference  in  average 
probability  of  detection  from  one  model  to  another  is  very-  noticeable.  Differences  in  model 
output  is  as  much  as  34 n<  predicting  the  performance  of  the  same  pattern.  Since  there 
was  not  much  difference  in  EPOD  between  the  five  patterns  using  one  model,  the  number  of 
patterns  considered  was  reduced  to  two.  Only  the  circle  and  4-4-4-4  pattern  were  pro¬ 
cessed  in  subsequent  scenarios. 
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Normal  Case 


l  ime  late  is  a  factor  for  I  POD  in  the  normal  case.  The  time  it  takes  an  aircraft  to 
arrive  on  station  changes  the  distribution  for  initial  target  location.  At  time  late  equal 
zero  TASDA,  AZoi,  and  sl’A.M  all  have  the  same  normal  target  location  distribution  as 
shown  in  Figure  11.  At  .‘10  min  lime  late,  the  normal  distribution  spreads  out  as  shown 
in  Figure  12  and  figure  12.  After  .‘10  min  the  probability  that  the  target  is  at  the  original 
estinuiled  position  is  smaller.  Since  the  target  location  distribution  does  change,  the 
eventual  probability  of  detection  values  are  different  for  each  time  late.  These  values  are 
given  in  Appendix  D  under  Normal  Case.  In  this  study  normal  case  output  CPOD  values 
were  averaged  over  three  Foils  and  three  frequencies  and  over  each  time  late.  Table  7 
Minus  averaged  values.  The  three  sets  of  output  show  that  AZOI  is  consistently  the  most 
optimistic  of  the  three  programs.  This  is  a  change  from  the  uniform  case  where  SPAM 
was  the  most  optimistic.  As  in  the  uniform  case,  the  averaged  values  were  very  close 
for  different  patterns  using  the  same  model.  The  ellipse  and  chevron  skew  patterns  were 
processed  subsequent  scenarios. 


Scenario  2 

Cumulative  probability  of  detection  was  calculated  for  two  patterns  with  a  uniform  and 
normal  initial  target  location  distribution.  The  first  flight  arrives  on  station  after  one 
hour  time  late  ;uul  remains  on  station  for  four  hours.  A  standard  deviation  in  target  speed 
O  introduced.  The  last  three  flights  arrive  on  station  with  no  time  late  and  remain  on  sta¬ 
tion  for  four  hours.  I  his  scenario  is  simulated  bv  a  computer  run  with  one  hour  time  late 
.uni  PI  hours  on  station. 

Tin*  values  in  I  able  s  were  obt .lined  b\  examining  CPOD  curves  and  attaching  a  nunneri¬ 
es  1  v  ,.lue  to  each  curve.  In  order  to  arrive  at  a  numerical  value  the  slope  of  the  curve  and 
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DISTANCE  (nm) 


Figure  12.  Normal  Target  Location  Distribution  for  TASDA 
and  AZOI  at  30  min  Time  Late. 


DISTANCE  (nm) 


Figure  13.  Normal  Target  Location  Distribution  for  SPAM 
at  30  min  Time  Late. 


its  maximum  value  should  he  considered.  A  numerical  value  based  on  area  under  a  curve 
was  used  since  area  under  a  curve  depends  on  rise  time  and  curve  height.  For  example,  in 
Fiitn re  11,  Curves  A  and  15  rise  to  (he  same  probability  of  detection  equal  to  20.  s.  However, 
curve  V  rises  to  that  I’nl)  faster  than  curve  I!  and  has  a  higher  numerical  indicator  attached. 
In  Figure  Id,  curve  C  rises  as  fast  as  curve  A,  but  curve  C  rises  to  a  smallr  r  I ’f  >f )  maxi¬ 
mum.  Therefore,  the  tactic  associated  with  curve  \  has  a  higher  measure  of  effectiveness 
than  curve  C  and  thus,  has  a  hit; her  numerical  indicator. 

The  area  under  the  (  Pol)  curves  was  estimated  b\  visually  inspecting  the  curves.  The 
method  used  to  calculate  the  CPOI)  curve  indicator  v  alues  is  shown  in  Figure  C-2  in  Appen- 
dix  C. 

I  niform  Case 

C!’(  >1)  indicator  values  for  the  uniform  case  were  averaged  over  three  ]  ( (Ms  and  three 
trequencics,  and  the  averages  are  presented  in  ruble  s.  lhe  indnator  values  lor  each 
curve  are  shown  in  Appendix  C.  In  this  scenario  onlv  two  patterns  were  considered  and 
computer  output  showed  only  partial  agreement.  However,  it  is  important  to  remember 
Unit  in  the  uniform  ease  initial  target  position  is  unknown.  For  all  three  models  indicator 
values  fell  within  a  few  percentage  values  of  each  other  for  the  1  — 1—4—1  and  circle  patterns. 
However,  the  difference  between  models  is  as  much  as  ;>7  in  the  1)-I)  case  and  i>  notice¬ 
able  in  all  cases. 

Normal  Case 

Cl’i  >1)  indicator  values  for  l lie  normal  case  were  averaged  over  three  FOMs  and  three 
frequencies  and  the  averages  are  shown  in  Table  8.  lhe  indicator  values  are  shown  in 
Apnvndi'-  c.  In  this  case,  initial  target  location  points  arc  clustered  toward  the  center  of 
the  ■  'per  lion  area  for  F  A-d >.\,  a /( >i  and  SIMM.  One  would  expect  these  three  models  to 
evaluate  li  e  na-'cins  imilarh  since  there  i<  now  some  control  on  target  location.  I  able  8 
- i •  * > v  av  ef-u.c  I'i'ob  il)> ! it v  ot  defection  values  remttiti  noticeably  different.  Note  that  the 
a  .Mtrr.a.  ■  u  Herein  e  !-  i  i  ><<i  ’in  S-S  deptii  with  the  ellipse  whereas  in  the  uniform 
>•  ■ :  f  '  mi. a.  fli  fl'ei  er.i «  'a-  2-1  '<i  me  circle  aw  die  S-S  depth.  This  relationship 

'  v  differences  can  be  observed  for  hot h  patterns  at  nil  depths.  The  long  on  station 
t  i  *  ■  '  F  '"ii'is  appear.-.  to  lure  an  a  fleet  on  t  lie  'urge  differences. 


'  "  '  •' . '  c  pi  i dvil hi  it  v  ->t  e<»! notion  casual  was  calculated  for  two  paiterns  with  a  uni- 
;  i  >.  iiiiil  i:  H  .:•!  t  ■  rget  location  dist  ri'mtior.  l  our  Piiglns  arrive  on  station  with  l 
•  ,-t  at  o  i - •  : 1 1 1 i t .  on  station  for  four  hours.  !  ach  nigh'  searches  :t  quadrant  of 

A  ■  ’  Held.  Siam  tin  sonobtiov  pattern-  a  »•<  a  mint -trie  with  respect  to  the  center 

a  ,  t  acli  .(uadrajit  has  an  hmniical  so.iobuov  geometry.  Therefore,  only  computer 

-in.  !  i,o  on  ont  .(uudrant  were  necessarv.  h.  stead  of  obtaining  C  Hnl)  indicator  values 
■  :  '  ••  -p.  i "tin v  cur  es,  av  •  rage  v  alues  were  obtained  bv  dividing  the  ( '  1  ’*  -I >  curves 

•  ■•'o  t  tu.e  -i  (■;*-•  '-shown  in  i'i  uii  ■  .  of  Appendix  1  .  ('  P(  > I  >  indicator  values  lot  the  uni- 
!"!>!•  ■:  *  *  n  I  case  w«  r<-  ,vi  r  wed  over  three  F<  M.-  and  three  frequencies  and  are 
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Figure  11.  .mes  A  and  B:  Numerical  Value  for 
Pattern  Selection  by  Rise  Time. 
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Figure  1.1.  Curves  A  and  C:  Numerical  Value  for 
pattern  Selection  by  FPOD. 
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Percent  Probabililv  of  Detection  for  Scenario  3 
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In  the  uniform  ease  results  in  Table  i)  show  that  there  are  still  significant  differences 
in  mode!  outputs  at  all  depth-..  However,  shortening  the  on  station  time  appears  to  have 
drought  the  values  closer  together  in  the  normal  case.  There  is  only  2.3'  maximum  dif¬ 
ference  at  the  S-s  depth. 


Scenario  I 

The  circle  and  chev  ron  skew  (Appendix  At  l»toy  fields  were  considered.  1  he  environ¬ 
ment  was  i  (tinged  from  the  Pacific  to  the  Atlantic  for  a  different  type  of  comparison.  I'he 
so’.ui.l  speed  profiles  in  I  igurcs  2  and  •'!  show  that  there  is  no  CZ  so  sound  propagation  i> 
bv  direct  path.  The  propagation  loss  curves  in  Appendix  B  indicate  bottom  bounce  is  present 
with  a  I  -inn  ft  .,<>mvo  receiver  depth  combination.  Values  in  Table  10  show  average  p<  d) 
'•tdues  ebiained  'or  the  AtlaiHie  cm  i  ronment .  llvdrophune  depths  of  100  and  1000  ft  were 
chosen  foi  simulation  and  the  results  art  shown  in  Table  10. 

fV  ,  sv cnaiie  wa-  .  liesen  t « »  show  how  th.e  models  respond  in  a  d.iect  path  environment. 
A-.  <  \f'  ••>(-d,  .(  .ah  os  a  re  I  >wcr  than  in  the  other  scenarios  because  of  the  absence  of  a 

convert  1  m  <■  /.one  on  inonmc:r.  the  one  noticeable  result  jn  the  output  averaged  in  Table 
10  i-  i  ho  hi  ah  ev  aluatmn  bv  AZ(  i|  using  the  chevron  skew  pattern  in  the  normal  case.  Actu¬ 
ally,  thi  .  i  «  '■  :  one  w  mid  prediv  t  for  i  he  chev  ron  skew  pattern  since  the  cov  erage  is 
greatest  it  the  center  where  the  target  is  most  likely,  i  he  other  two  models,  however, 
did  not  ‘•.crease  fiieii  prediction  -  -  igni  ficam  I\  from  I  he  uniform  case. 


Table  10 

Percent  1’rolabilitv  of  Detection  for  Scenario  A 


CNIKOKM 

CASK 

s-s 

D-D 

01  ITS 

1  II  KAM 

ClliCl  1 

(TASDA  1 

id.  i 

s.  7 

CHIC  K1 

t  \/.oi> 

2S.  () 

2S.  i 

!  CIRC  LI 

(SPAMt 

2s.  7 

2 5.  7 

SKI  >PE 

\v  \t  kr 

Clb’C  I  I- 

(1.  0 

s.-l 

cult  i,i 

21 . 1! 

23.  7 

CIRCI.1 

:’.g.  1 

30.  7 

NORMAL  CASK 

S-S 

D-D 

|  CiCI.K  STR  KA.tf 

CHEVRON  SKEW 

2s.  s 

12.3 

CHEVRON  SKEW 

67.2 

30.  0 

CHEVRON  SKEW 

2-1.  1 

32.  0 

SLOPE  WATER 

CHEVRON  SKEW 

10.  1 

0.  1 

CHEVRON  SKEW 

41.  s 

3(1.  1 

CHEVRON  SKEW 

22.  0 

32.  .1 

VI.  COXC  I.KSloNS 

I  \peeted  probn!>ilitv  of  detection  is  calculated  using  four  different  methods.  ADI  l’S 
calculates  the  average  conditional  probability  in  a  symmetrical  subsection  of  a  buoy  field 
pattern  where  the  submarine  is  assumed  to  be  a  random  target.  TASDA  computes  prob- 
abilin  of  detection  br  one  or  more  b uovs  using  Monte  Carlo  simulation  of  100  random 
tracks  through  a  buov  field  and  a  detect  (POD  ;  1)  or  no  detect  (POD  -  Of  analysis  for 
each  track.  A/< )!  uses  a  binominal  probability  function  to  calculate  KPOD  where  a  value 
for  time  along  a  track  is  chosen  randomly.  SPAM  uses  a  Monte  Carlo  simulation  and 
calculates  I  pi  >D  bv  averaging  a  probability  obtained  through  assuming  independence  of 
bum  s  with  a  probnbilitv  olnained  by  assuming  complete  dependence  of  buoys,  since  the 
method-  of  calculation  are  different  in  each  program,  the  output  for  KPOD  is  not  exactly 
the  same. 

ADI  l’s,  TASDA,  AX<  d  and  SPAM  all  calculate  cumulative  probability  of  detection  as 
a  function  of  time.  However,  l he  method  of  calculation  is  different  in  all  four  programs. 
ADI  PS  uses  an  exponential  expression  to  estimate  a  conditional  cumulative  probability 
pit  i  as  discussed  in  Section  l.  TASDA  uses  a  detect  -  no  detect  Monte  Carlo  methodology. 
AZ(  >1  uses  I  Pi  >1)  calculation  to  obtain  tin  estimate  for  the  lower  limit  of  CPOD.  In  the 
case  where  the  target  location  is  assumed  to  have  a  uniform  distribution,  the  C POD  does 
not  change  with  time  and  i.-,  equal  to  I  p<  iD.  In  the  normal  case,  KPOD  and  CPOD  are  dif— 
IV* rent  since  calculations  for  the  normal  case  are  time  dependent.  SPAM  uses  a  cumulative 
time  window  to  limit  Pi  >1)  calculations  for  CPOD.  Since  different  techniques  are  used  for 
each  program,  it  is  not  expected  that  CPOD  output  will  be  exactly  the  same  for  each  pro¬ 
gram.  ADI  I’s  outputs  P  and  P( t > ;  TASDA  outputs  CPOD,  KPOD  and  mean  time  to  first 
detection  (MTKDi:  AZ<  >!  outputs  CPOD,  1  Pol),  standard  deviation  (SIC.MAi  and  probnbil¬ 
itv  that  a  target  will  be  detected  given  that  it  has  not  been  detected  in  (lie  past  (P  I ) I  'NDTt; 
spAM  outputs  Cpol)  and  I  pop.  Although  the  svmbols  CPOD  and  KPOD  in  Table  ii  are 
used  for  TASDA,  A/.<  d,  and  SPAM,  the  output  for  these  parameters  may  not  be  the  same. 


m  K d li .wiim  ci 'iK'lus ii ms  art  indicated  I)'.  this  1 1 n .* \ . 


!,  TASDA  i-  ( onserv  :ili\ e  at  l.mv  l'(  >\Is 

Inlike  A  X »  1 1  or  s|’AM,  TASDA  vf<  iui  rrs  contact  during  two  consecutive  time  steps 
before  a  lit  Uction  is  made.  TASDA  const  met  s  prolialiilit\  areas  around  each  buov  ijounded 
In  tv:. .  i-oiu'intric  circles  as  in  ligure  K>.  I  (u  se  areas  represent  the  locations  where  the 
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1  inure  1C.  Target  I'raok  for  IASI )  A  and  SI  >  u  I  teuton  Around  a  Sonobuov.  I>(  il)  =  1. 

signal  excess  St  is  greater  than  zero.  It  the  target  remains  in  one  of  those  areas  during 
two  consecutive  time  steps,  then  the  l><  i|)  1  tor  that  time  period.  <  itherwise,  the  l)(iD  = 

There  will  ho  a  set  of  these  values  for  each  of  the  100  '|  ASDA  trials.  The  values  are 
simnit'd  cv  ei  100  trials  to  conie  up  with  ('!'<  i|)  values  for  each  time  interval.  At  low  TttMs 
the  areas  hounded  bv  concentric  circles  become  ver\  smtil!  and  as  a  result,  the  C'l’ni) 
values  drop  rapidly  to  zero.  A7.(  »I  and  SPAM  compute  P<  >P  versus  range  for  each  buoy 
rather  than  assigning  a  0  or  1  value  tit  a  particular  ramie. 

1’.  SPAM  is  Optimistic  at  low  i’l'Als 

If  the  taruet  's  at  a  rariuc  >uch  that  a  faint  line  will  be  created  on  a  uram,  then 
there  trust  be  some  .amount  of  time  during  which  the  taruet  remains  at  a  detectable  range 
in  order  for  the  operator  to  me  a  line;  call  it  t0.  Keen  aifer  tQ  time  has  passed,  the  oper¬ 
ator  nav  no<  see  the  line.  The  prnbabjlitv  that  the  operator  sees  the  line  increases  up 
until  -o  we  time,  t7.  Iho'-e  two  values  for  time  tire  incorporated  into  a  cumulative  time 
window  t('TW  i  model  which  examines  the  probabilities  between  these  times  in  increments 
of  time  i],  I  lie  smaller  the  value  of  ty,  the  higher  'he  POD  value.  Standard  operating 
values  fot  <PAM  are  Up,  tj,  <z,  -  f. 2.1,  . 50|.  This  means  that  the  Target  must  re¬ 

main  at  a  detectable  ranee  f<w  onlv  Id  minutes  If)  achieve  the  probnbil  tv  o'  detection  pro  - 
dieted  b\  tie  standard  passive  sonar  equation.  The  v  alues  to,  1  j ,  p-  in  the  (’T\V  caused 
ptwdielio  ■>.  for  SPAM  to  be  optimistic  at  low  l'oMs.  Detection  can  occur  at  any  time 
before  or  a  tier  Ip,  depending  on  chance. 

.'t.  SPAM  Results  are  more  Conservative  at  High  1'uMs 

I'iie  ( '  r\V  model  builds  detection  probal'ilities  experimentally  for  a  target  re- 
mainitu  at  a  fixed  rang*-  uni ; !  tins*  ’7  with  the  probability  predicted  by  the  sonar  equation 
achic  •  t  .titf.r  time  t„  has  pa  ;s«!.  This  function  a. so  smoothes  the  jump  from  low  prob- 
tibilitie  associated  with  iow  liO.is  to  hi  :h  prob.d  Eli;  ie«  associated  with  high  TTiMs. 


I.  I  \SDA  and  AZ<  )I  Show  a  lump  in  (’POP  as  l-'OM  Increases 


A  normal  distribution  around  the  figure  of  merit  is  used  in  the  SPAM  model. 

The  belief  is  that  the  area  under  a  normal  curve  is  equal  to  the  probability  of  detection. 

SI  =  0  when  FuM  -  PI..  At  that  point  there  is  a  50r.'  POP.  As  shown  in  Figure  17, 
spAM  POD  gradually  increases  ns  SF  increases.  Truncation  has  an  effect  on  CPOI)  values 
for  AZof  with  low  FOMs.  For  low  FOMs  there  is  very  little  signal  excess  and  probnbilitv 
depended  on  the  left  half  of  the  curves  above;  therefore,  the  CPOD  values  dropped  quicklv 
to  0.  TASPA  uses  a  Hician  Distribution  to  simulate  short-term  fluctuation  in  the  PI.  curve 
and  calculations  are  based  on  the  assumption  that  below  the  FOM  an  operator  is  unlikely 
to  detect  a  target.  Detection  is  either  YFS  or  NO  based  on  SF  0  after  employing  the  Tlician 
fluctuations.  As  is  apparent  from  the  graphs  in  Figures  17,  18  and  19.  AZOI  and  TASDA 
have  probability  distributions  which  drop  off  quickly  to  zero  below  the  FOM.  Therefore, 
TASDA  and  AZOI  exhibit  breaking  points  when  the  FOM  reaches  a  small  enough  value.  At 
these  FOMs  (  POD  falls  rapidly  to  zero. 


r>.  AZOI  Shows  Small  Variance  Between  Patterns 

AZOI  differs  from  TASDA  and  SPAM  in  that  TASDA  and  SPAM  use  successive  time 
steps  along  a  target  track  for  a  particular  trial.  After  a  track  has  been  generated  for  a 
particular  trial,  AZOI  picks  random  points  along  that  track  and  evaluates  POD  values  at 
those  points.  Therefore,  the  time  increments  are  not  evenly  spaced  but  randomly  spaced. 

At  each  (joint,  the  range  to  each  buoy  is  calculated  and  the  pi’obabilities  associated  with 
each  buoy  are  summed  to  give  a  POD  for  the  target  at  a  point  along  the  track.  In  this 
manner,  all  contributions  arc  considered.  The  probability  at  a  point  along  the  track  is  not 
zero  or  one  as  with  TASDA,  nor  is  execution  terminated  when  a  detection  is  encountered  as 
with  SPAM.  So  this  method  yields  more  values  for  POD  along  a  track;  and  a  cumulative 
time  window  is  not  employed  as  is  the  case  with  TASDA  and  SPAM. 

These  five  differences  in  the  three  computer  simulation  models  result  from  various 
interpretations  of  the  tactical  sonobuoy  situation  under  consideration.  These  interpretations 
result  in  different  output  values  for  each  of  the  three  simulation  programs. 

The  large  differences  between  computer  models  require  further  investigation.  Future 
studies  should  be  based  on  a  real  world  application  of  the  sonobuoy  programs.  Heal  world 
mission  planners  usually  have  SOS  IS  or  other  intelligence  information  for  input  to  the 
model  and  are  interested  in  POD  versus  buoy  spacing.  The  output  data  sets  in  this  report 
are  not  designed  to  reflect  a  real  world  application  of  the  sonobuoy  programs,  but  acquaint 
the  reader  with  the  amount  and  type  of  data  output  generated  by  considering  a  few  pattern 
configurations  under  several  different  conditions. 
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SIGNAL  EXCESS  (dB| 


Figure  17.  SPAM:  POP  vs  SK  Used  in  the  Calculation  of  a  Lateral 
Range  Curve  (Normal  Distribution). 
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Figure  18.  AZOI:  POD  vs  SL  Used  in  t lie  Calculation  oi  a  Lateral 
Lange  Curve  (ROC  Curve). 


?4  70  16  I?  8  4  0  4  8  1 7  16  .'0  ?4 

SIGNAL  EXCESS  IdBI 


Figure  1(>.  i  \  'DA:  POD  vs  SI  I 'sc I  in  the  Calculation  of  a  l  ateral 
•  Lumet  uive  (Detect-No  Detect  Assumption). 
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Figure  A-l.  Circle  pattern  used  in  scenarios  1,  2,  3,  and  -4 
for  a  uniform  target  distribution. 


Figure  A ~2.  Ellipse  pattern  used  in  scenarios  1,  2,  3,  and  4 
for  a  normal  target  distribution. 
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f  igure  A-3.  4-4-4-4  pattern  used  in  scenario  1  for  a  uniform  and 
normal  target  location  distribution  and  in  scenarios  2, 
3  and  4  for  a  uniform  distribution. 
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Figure  A-4.  chevron  or  X  pattern  used  in  scenario  1 
for  a  uniform  target  distribution. 


figure  A->.  chevron  skew  pattern  used  in  scenarios  1,  2,  2,  and  4 
for  a  normal  target  distribution. 
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figure  A-G.  5-G-3  pattern  used  in  scenario  1  for  a  uniform 
and  normal  target  location  distribution. 
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Figure  A-7.  Brushtac  pattern  used  in  scenario  1  for  a  uni¬ 
form  target  distribution  and  for  a  normal 
target  distribution. 
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Figure  13-2.  Propagation  loss  curves  at  a  sh;illo\v  receiver  depth  and  deep 
source  depth  for  scenarios  1,  2,  and  3. 
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B-4.  Propagation  loss  curves  exhibiting  direct  path  sound  transmission 
with  no  convergence  zone  used  as  input  to  the  tactical  models  in 
scenario  4. 
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Thi>  appendix  contains  n  summarv  of  the  computer  runtime  for  each  ol  the  scenarios. 

[t  should  I'e  no!ed  that  the  runtime  for  AX<  >1  includes  the  calculation  ol  P  DP  NDT  in  addi- 
tjonto  I  PtiDttnd  t'Pnl).  Also,  AZ(  i|  includes  pattern  processing  time.  Runtime  must 
! >e  miltiplied  In  a  factor  of  <10  when  estimating  minicomputer  runtime.  Therefore,  the 
values  can  t>e  considered  in  hours  rather  than  minutes  when  making  minicomputer  run¬ 
time  estimates. 

Also  in  tliis  appendix  are  the  original  data  generated  bv  the  programs  f  Figures  C'-l 
through  The  data  are  presented  in  graphical  form.  The  original  data  were  reduced 

to  single  I  POD  and  ('POD  indicator  values  for  the  purposes  of  program  comparison  in  this 
report.  The  method  for  obtaining  the  indicator  values  is  shown  in  Figures  C-ti.'J,  C'-fil, 
('-<)■).  The  derived  indicator  values  for  eacli  POM,  frequency  and  depth  combination  are 
listed  in  the  tallies  in  the  appendix.  These  values  were  averaged  to  obtain  values  for 
analysis  in  the  text  of  this  report. 
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Table  (  -1.  Runtime  Psinn 
The  I'nivae  ll()s  Computer 


(  IT  Time 

l/( )  Time 

Total  Time 

Scenario  1 

TASDA 

til.  7  min 

1.  ">  min 

2  min 

,\  yx  >i 

1  si).  ;i  min 

.  tj  min 

iso.  i,  min  i 

SPAM 

lit).  0  min 

1.  j  min 

1 1 1 .  j  min 

x 

Scenario  2 

I  ASI  'A 

17.  t)  m  in 

1.  2  min 

1.2  min  1 

A /i »! 

hi.  0  min 

.  ti  min 

iti.  :>  min 

SIAM 

2t>.  ti  min 

0.  !)  1’ i it! 

21.  2  min 

| 

Sivnnrio 

1 

I'ASI  >.\ 

U .  (t  min 

1. 1!  min 

7.  !*  rain 

A/.i  d 

H .  u  min 

.  ti  min 

M.  ti  min 

Sl’A.M 

12.1  min 

.  St  min 

lti.  0  min 

S<  ■•.'ll*  (fit;  \ 

. 

I’ASMA 

min 

2  min 

7.  0  min 

:  A/.OI 

1  2.  7  min 

.  ti  min 

12.  0  min 

SIAM 

ll.fi  mi.n 

.  7  min 

12.  .7  min 

_ _ _ i 

''lie  output  u  -ed  10  const  ruct  graphs  of  the  original  data  w  as  obtained  from  compute' 
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I'i^ui'e  C  15.  Comparison  of  TASDA .  A/.OI  ;in<l  SPAM  with  a  uniform  initial 
target  location  at  30  minutes  time  late  with  a  deep  tarnot  and 
a  deep  receiver  in  seenario  i  usinii  the  chevron  pattern. 
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SCENARIO  4  -  UNIFORM  CASE 
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Figure  C  -  57 .  Comparison  of  TASDA ,  AZOI  and  SPAM,  with  a  uniform  initial 
target  location  at  one  hour  time  late  with  a  receiver  frequency 
of  100  Hz  in  the  slope  water  environment  in  scenario  -1  using 
the  circle  pattern. 
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Figure  C  58.  Comparison  of  l’ASDA .  A  7.01  and  SPAM ,  with  a  uniform  initial 
target  location  at  one  hour  time  late  with  a  receiver  frequency 
of  >100  Hz  in  the  slope  water  environment  in  scenario  1  using 
the  circle  pattern. 
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Figure  C-59.  Comparison  of  TASDA,  AZOI ,  and  SPAM  with  a  normal  initial 
target  location  at  one  hour  time  late  with  a  receiver  fre<|ueney 
of  100  Hz  in  the  Gulf  Stream  environment  in  scenario  -1  using 
the  chevron  pattern. 
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Figure  C-60.  Comparison  of  TASDA.  AZOI.  in  !  SPAM  with 
target  location  at  one  hour  time  late  with  a  r< 
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Figure  C-63.  Scenario  1  -  Kventual  Probability  of  Detection  for  a 

shallow  target  and  receiver.  The  eventual  HOD  value 
is  t!  ■*  maximum  l'i  HOD  attained  by  the  curve.  Kventual 
HOD  is  obtained  by  observing  the  maximum  HOD  for 
each  CHOD  curve.  These  maximum  values  are  summed 
over  frequency  and  FOM  to  give  a  percent  HOD  for  a 
uniformly  distributed  shallow  target  and  a  shallow 
receiver.  Not  for  example:  (.14  +  .50  +  .96  -*  .  22  + 

.58  +  .86  4  .21  + .26  •  ,76)/9  =  49.9. 
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Table  t'-3:  Scenario  1 
Kventual  Probability  of  Detection 
for  i  Uniform  Target 

at  30  min  Time  I.ate:  Shallow  -  Deep  Case 
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I  aide  :  Scenario  1 
I  ventnal  I'rolj.ibiliU  of  Detection 
for  a  Normal  Target 
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Table  C-9:  Scenario  1 
Eventual  Probability  of  Detection 
for  a  Normal  Target 

at  1  Hr  Time  l.ate:  Shallow  -  Deep  Case 
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Table  (.'-10:  Scenario  1 
Kventual  Probability  of  Detection 
for  a  Normal  Target 
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Table  (  -12:  Scenario  1 
Lventiuil  Probability  of  Detection 
for  a  Normal  Target 

at  4  hrs  lime  Late:  Shallow  -  Deep  Case 
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Table  C-12:  Scenario  1 
Eventual  Probability  of  Detection 
for  a  Normal  Target 
at  4  Hrs  Time  Late:  Deep  -  Deep  Case 
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Figure  C-0-1.  Scenario  2:  Average  Probability  of  Detection  a  shallow 
target  ami  receiver.  The  average  POD  value  is  obtained 
b\  visually  estimating  the  area  under  the  C'PoD  curve. 
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Table  C-13:  scenario  2 
Average  Probability  of  Detection 
for  a  Uniform  Target 
at  1  Hr  Time  Late 
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Average  Probability  of  Detection 
for  a  Normal  Target 
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Figure  065.  Scenario  3  and  4:  Average  Probability  of  Detection  for  a 
shallow  target  and  receiver.  The  average  POD  value  is 
obtained  by  summing  the  POD  values  at  each  time  incre¬ 
ment  .and  dividing  by  the  total  number  of  time  increments 
to  estimate  the  area  under  the  CPOD  curve.  The  average 
values  are  summed  over  frequency  and  FOM  to  give  a 
percent  POD  for  a  uniformly  distributed  shallow  target 
and  a  shallow  receiver. 
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Table  ('-  In :  Scenario 
Average  Probability  of  Detection 
for  a  Normal  Target 
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Table  C-l*:  Scenario  4 
Average  Probability  of  Detection 
for  a  Normal  Target 
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